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Abstract: We compute the two- loop enhancement factors for our earlier one-loop 
calculations of leading {1/Q) power corrections to the mean values of some event 
shape variables in deep inelastic lepton scattering. The enhancement is found to 
^| be equal to the universal "Milan factor" for those shape variables considered, pro- 

vided the one-loop calculation is performed in a particular way. As a result, the 
. ^ ' phenomenology of power corrections to DIS event shapes remains largely unaffected. 
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1. Introduction 

In a recent paper [1] we presented a theoretical study of the leading power-suppressed 
corrections to the mean values of various event-shape measures in deep inelastic scat- 
tering (DIS). Our investigation was based on an analysis of one- loop Feynman graphs 
containing a massive gluon, which probes the sensitivity of different observables to 
long-distance physics [2, 3, 4] and is closely related to the analysis of infrared renor- 
malons [5, 6, 7, 8]. We found that, as in e + e~ annihilation, the leading corrections to 
DIS event shapes are normally of order 1/Q. This prediction has since been verified 
experimentally [9]. 

By making the further assumption of universal low-energy behaviour of the 
strong coupling, we also estimated the magnitudes of 1/Q corrections to event shapes 
in DIS relative to those observed in e + e" annihilation [10, 11]. Here the one-loop 
massive-gluon analysis suffers from some ambiguities and deficiencies. The ambigu- 
ities arise from the sensitivity of some event shapes to the way in which the gluon 
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mass is included in the definition of the shape variable, the phase space and the 
matrix elements. This problem is related to the inadequacy of the massive-gluon 
approach pointed out by Nason and Seymour [12]: unlike the total e + e~ annihilation 
cross section, for example, event shapes are sensitive to the precise way in which 
a virtual gluon fragments into observable particles. In other words, they are not 
sufficiently inclusive to be expressible in terms of the distribution of an effectively 
massive gluon. 

In a two-loop analysis of this problem in e + e~ event shapes, however, the Mi- 
lan group [13] has found that, for a wide class of shape variables, the effects of 
non-inclusiveness amount simply to an enhancement of the "naive massive-gluon" 
estimate of the 1/Q correction by a universal "Milan factor", with numerical value 
M. ~ 1.49. The two-loop analysis also clarified the way in which the gluon mass 
should be included at the "naive" level, in order that universality should be mani- 
fest. 

In the present paper we perform the same type of two-loop analysis for event 
shape variables in DIS. We find a universal enhancement of the 1/Q corrections to 
the current jet thrust, mass and C-parameter, with the same enhancement factor 
Ai, provided the "naive massive-gluon" estimate is computed in a particular way. 
This resolves the ambiguities we encountered in the one-loop calculation of power 
corrections to certain DIS shape variables. 

In the case of the current jet broadening variable, which we also studied at one- 
loop level in Ref . [1] , there are additional complications arising from recoil contribu- 
tions [14], and therefore we postpone the two-loop analysis to a future publication. 

The layout of the present paper is as follows. In Sect. 2 we define convenient 
kinematic variables for the study of multi-parton emission in DIS. These variables 
are used in Sect. 3 to express the various event shape variables that we investigate. 
We shall argue that the dominant two-loop contributions to the leading power cor- 
rection come from the kinematic region in which the emitted gluons are soft. The 
contributions from this region are described in Sect. 4. Sect. 5 explains the 'disper- 
sive' method we use to calculate the power corrections. The actual calculations are 
then performed in Sect. 6, and finally the results are summarized and discussed in 
Sect. 7. 



2. Kinematics 

It is convenient to write the momenta k{ of radiated partons (gluons and/or quark- 
antiquark pairs) in terms of Sudakov (light-cone) variables as (Fig. 1) 

h = a t P + PiP' + k tl , OiPi = k 2 JQ 2 (2.1) 
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Figure 1: Deep inelastic scattering with multi-parton emission. 



where P and P' are light-like vectors along the incoming parton and current direc- 
tions, respectively, in the Breit frame of reference: 

P = xP p , P' = xP p + q, 2(P ■ P') = -q 2 = Q 2 (2.2) 

where P p is the incoming proton momentum and x = Q 2 /2(P p ■ q) is the Bjorken 
variable. Thus in the Breit frame we can write P = \Q{1, 0,0,-1) and P' = 
\Q(l, 0, 0, 1), taking the current direction as the z-axis. Particles in the current 
hemisphere H c have > a« while those in the proton remnant hemisphere H r have 
on > fa. 

The momenta p and p' of the initial- and final-state quarks can also be resolved 
along the Sudakov vectors P and P'. From momentum conservation, 

p = (1 + a' + J2 a i) p 

p' = a 'P + (l-'£p i )P'-'£k ti (2.3) 

where 

«'(i-EA) = (E4) 2 /Q 2 - (2-4) 

The initial quark is assumed to be collinear with the proton direction and is therefore 
aligned along P. It may have an "intrinsic" transverse momentum of the order of 
the hadronic mass scale A ~ 1 GeV, but this will lead to corrections of order \ 2 /Q 2 , 
which we neglect. The final quark momentum p' is mostly along P'\ its component 
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a' along P is small. In perturbation theory, a' is of order a s (Q 2 ) on the average, 
while in the soft region (all ku ~ A) it is of order \ 2 /Q 2 . When some emissions 
are hard and some soft, the non-perturbative contributions to event shapes are at 
most of order a s (Q 2 )\/Q, which we also neglect. Thus, as far as X/Q corrections 
are concerned, a' can be regarded as contributing only to the perturbative part of 
event shapes, and we can compute the non-perturbative part from a study of the 
soft region alone. 

In the soft region, the matrix elements for multi-parton emission in DIS will be 
the same as those for e + e~ annihilation into quark- ant iquark plus soft partons at 
cm. energy-squared Q 2 , with the outgoing antiquark replaced by an incoming quark 
of momentum p. This means that we shall be able to take the results for various 
matrix element integrals directly from the study of e + e~ event shapes [13], adjusted 
according to the definition of DIS event shapes. 

3. DIS event shape variables 
3.1 Current jet thrust 

Consider first the current jet thrust Tq, defined by the sum of longitudinal momenta 
in the current hemisphere normalized to Q/2: 

We have assumed that the outgoing quark momentum p' lies in the current hemi- 
sphere, because the probability that this is not true is of order a s (Q 2 ), and the 
associated non-perturbative contribution will be at most of order a s (Q 2 )X/Q, which 
we neglect. In the notation of the previous section, we have 

r = 1 -T = EA-E(A -«*) + «' 

= E A + E a i + a ' = E min{«i, A} + ol . (3.2) 

Where not explicitly stated, the sums extend over all partons, irrespective of which 
hemisphere they are in. Here the contribution of those partons that are not directly 
in the current hemisphere represents the loss of momentum of the outgoing quark. 

Note that the expression (3.2) is additive, in the following sense: the mean 
value of the event shape receives a set of perturbative contributions, from a' and 
mm{ai, fa} for hard parton emissions, plus a non-perturbative contribution from 
soft emissions, which can be computed from the expression 

tq ~ J2 min i a i^i} » ( 3 - 3 ) 
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valid in the soft region, and added to the perturbative part. 

If instead we measure the thrust with respect to an axis chosen to maximize its 
value in the current hemisphere, we obtain the same expression in the soft region, 
since the transverse contribution is of second order in the small quantities ku/Q. 

Another possible way of defining the current jet thrust is to normalize it to the 
total energy in the current hemisphere, instead of Q/2. This gives 

T E = T Q /(1 - s) (3.4) 

where 

e = 1 - ^ £ E a ■ (3.5) 

The quantity e, which we call the energy deficit in the current hemisphere, is itself 
an interesting shape variable, not measured up to now. It is given by 

e = EA- £(A + ai)-af = ]TA-E^-«' • (3.6) 

H c H r H c 

and hence in the soft region 

r E = l-T E ^2Y,oc l . (3.7) 

3.2 Current jet mass 

The current jet mass is defined by 

P=(lLPa) IQ 2 ■ (3-8) 
\aeH c J 

In the notation of Sect. 2, 

p = L' + ^q,) . ll-^Aj -5>iA- (3.9) 

Once again we should ignore a' and terms of second order in ccj and/or fa when com- 
puting 1/Q corrections, in which case we have simply the additive soft contribution 

P ^ ( 3 - 10 ) 

3.3 C-parameter 

The C-parameter is 

C = 6 E \Pa\\Pb\sm 2 9 ab /Q 2 . (3.11) 

a,beH c 
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Note the factor of 6 replacing the 3/2 that appears in the definition of the C- 
parameter for e + e~ annihilation, because of the normalization to (Q/2) 2 in place 
of Q 2 . 

In the notation of Sect. 2, the general expression for the C-parameter is compli- 
cated, but in the soft region we obtain simply 

^-12£-^, (3.12) 

which is again additive. 

4. Soft two-loop contributions 

Consider now the contribution to a shape variable V when two soft, not necessarily 
collinear partons are emitted. Following [13], we write the two-parton phase space 
as 

da x da 2 d 2 k t \d 2 k t2 dad 2 k t d(j) 2 
ar 2 (&;i, k 2 ) — = —dzdm (4.1) 

OL\ a 2 7T 71 Ct 71 27T 



where 



a = a\ + a 2 , Oil — za , a 2 — (1 — .z)ck , 
m 2 = (fci + /c 2 ) 2 , fct = fcti + fc t2 (4.2) 

and is an azimuthal angle. The contribution to the mean value of V from one or 
two soft partons, which will be denoted by V, may be written in the form 

da d 2 k t ( a s (0) 



+ 2C F J dr 2 (h, k 2 ) Q|) 2 {C F M F + C A M A + n f M f )v 2 (h, k 2 ) , (4.3) 

where v±(k) is the contribution to the shape variable from a (massless) parton of 
momentum k and v 2 (k±, k 2 ) is the corresponding two-parton contribution. 

The first term on the right-hand side of Eq. (4.3) represents the lowest-order 
contribution of a single soft gluon (involving the ill-defined quantity a s (0)) plus the 
two-loop virtual corrections to it. The second term represents the contribution from 
the emission of two soft partons, Mp etc. being the parts of the relevant matrix 
elements-squared with the corresponding colour /flavour factors. Both terms are di- 
vergent, but we can rewrite V as sum of collinear and infrared finite parts by the 
following procedure. We introduce a term representing a contribution from the com- 
bined parton momenta, V\{k\ + k 2 ), and split the two-parton contribution into three 
parts: 

v 2 (h, k 2 ) = v^h+k^+lv^ki) + ui(fc 2 ) - ui(fci + k 2 ) ]+[u 2 (fci, h) - Ui(fci) - Mh) ] 

(4.4) 
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The first term defines the "naive massive-gluon" contribution, which treats the parton 
emission inclusively. The first expression in square brackets is the "non-inclusive" 
correction [13]. The second square bracket defines what we shall call the "non- 
additive" contribution. 

The combined momentum ki + k 2 is massive and therefore the definition of the 
contribution Vi(ki + k 2 ) is somewhat ambiguous. It does not matter how we define 
this quantity, as long as the expression used becomes equal to V\(ki) + ^1(^2) i n the 
soft and/or collinear (i.e. massless) limit. We define it as follows: if 

vi(ki) = v(ati, Pi) , a i p i = k 2 i /Q 2 , (4.5) 

then 

vxfa + fa) = v(a,P) , a = a! + a 2 , (3 = P1 + P2; &P = (k 2 + m 2 )/Q 2 . (4.6) 

The expression (4.3) may now be written as the sum of four finite parts, as 
follows [13]: 

V = V + V in + V ni + V na . (4.7) 

The first term is the so-called "naive" contribution. It incorporates only those two- 
loop terms proportional to the beta-function coefficient Pq = (11CU — 2n/)/3, which 
together with a s (0) build up the running coupling, as will be explained in the fol- 
lowing section: 

where f2 , the "naive trigger function", is given by 

ft (e) = / —v{a,P = e/a) . (4.9) 
Je a 

The second term in Eq. (4.7) is the "inclusive" correction, which combines the 
remaining non-Abelian part of the virtual correction with the corresponding "naive 
massive gluon" part of the two-parton contribution. This was shown in Ref. [13] to 
be given by 



V,„ = SC F C A [ **f J ^)\n 
J m 2 (kf + m?) \AtiJ 

■ n ((k 2 + m 2 )/Q 2 ) - n (k 2 /Q 2 )} (4.10) 

The "non-inclusive" correction is the additive part of the remaining non-Abelian 
(and quark-antiquark) two-parton contribution: 

V™ = ^ J dm 2 dk 2 t dzdct){^j 2 {C A M A + n f M f ){l ni {m 2 /Q 2 ,k 2 /Q 2 ,z,<j)) (4.11) 



k 2 (k 2 + m 2 ) 



m 4 
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where the "non-inclusive trigger function" is 

da 

l{k 2 +m 2 )/Q 2 a 



n ni (m 2 /Q 2 ,k 2 /Q 2 ,z,<j>) = I —[v 1 (k 1 )+v 1 (k 2 )-v 1 (k 1 + k 2 )}. (4.12) 

J(k 2 +m 2 )/Q 2 OL 



Note that the Abelian (Cp) part of this contribution is cancelled by virtual corrections 
[13]- 

Finally the "non-additive" correction is the rest of the two-parton contribution: 

V na = ^J dm 2 dk 2 dzd<f ) (^) {C F M F + C A M A + n f M f )Q na {m 2 /Q 2 ,k 2 /Q 2 ,z,<f)) 

(4.13) 

where the "non-additive trigger function" is 

Sl na (m 2 /Q 2 ,k 2 t /Q 2 ,zA) = l\ —[v 2 (k 1 ,k 2 )-vi(k 1 )-v 1 (k 2 )]. (4.14) 

J (k 2 +m 2 )/Q 2 OL 

All the quantities considered in Sect. 3 are additive, in the sense that v 2 (ki,k 2 ) = 
vi(ki) + vi(k 2 ) in the soft region, and so there is no contribution V na in these cases. 

5. Dispersive calculation of power corrections 

Our method for estimating power corrections will assume that the QCD coupling 
a s (k 2 ) can be defined down to arbitrarily low values of the scale k 2 and that it 
has reasonable analytic properties, i.e. no singularities other than a cut along the 
negative real axis. It follows that one can write the formal dispersion relation 

*s(k 2 ) = - r p s {m 2 ) (5.1) 



lo m 2 + k 2 

where the spectral function p s represents the discontinuity across the cut, 

Ps(m 2 ) = -i-Disc{a s (-m 2 )} = ^- {a s (m 2 e™) - a s (mV^ )} . (5.2) 

Our strategy will be to rewrite the results of the previous section in terms of 
p s (m 2 ) and then to study the effect of a non-perturbative modification to this function 
in the soft region. We eliminate a s (0) from Eq. (4.8) using Eq. (5.1), and insert 
p s (m 2 ) in the second-order terms in place of a 2 using the result 

The appropriate scale for the argument of p s in this substitution is m 2 , since the two- 
loop contribution in Eq. (4.8) generates the running coupling. We assume that the 
same argument can be used in Eqs. (4.10) and (4.11). Then, for example, Eq. (4.8) 
becomes 

V = f^p s (m 2 )Fo(m 2 /Q 2 ) (5.4) 
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where the naive characteristic function To{m 2 / 'Q 2 ) is 



Fo{m 2 /Q 2 ) 



7T JO 



n ((k 2 + m 2 )/Q 2 ) Q (k 2 /Q 2 ) 



k? + m 2 



k? 



(5.5) 



Similarly we may define Ti n etc. for the inclusive, non-inclusive and non-additive 
corrections, 



^n(m 2 /Q 2 ) 
T ni {m 2 /Q 2 ) 
T na [m 2 IQ 2 ) 
so that 



2C fC a 



Cpm 



alk 2 



kf + m 2 



In 



k 2 (kf + m 2 ) 



rrr 



{Q ((k 2 + m 2 )/Q 2 )-Q (k 2 /Q 2 )} 



4tt 2 A 
C F m 2 



r f dk 2 dzd<P (C A M A + n f M f ) Q ni (m 2 /Q 2 , k 2 /Q 2 , z, 0) (5.6) 
o J 



= I dk 2 t dzdct ) {C F M F + C A M A + n f M f )a na {rn 2 /Q\k 2 t /Q 2 ,z, ( t>)] , 



V = [^ Ps (m 2 )F(m 2 /Q 2 ) 
J m 2 



where T is the two-loop corrected characteristic function, 

3~ 3~ 3 ~t~ 3~ in ~t~ J~\ ni ~t~ 3~ na • 



(5.7) 



(5.J 



Non-perturbative effects at long distances are expected to give rise to a modi- 
fication in the strong coupling at low scales, 5a s , which generates a corresponding 
modification in the spectral function via Eq. (5.2): 



Sp s (m 2 ) = 7^-Disc |5a s (-m 2 )| 



(5.9) 



Inserting this in Eq. (5.7) and rotating the integration contour separately in the 
two terms of the discontinuity, we obtain the following non-perturbative contribution 
to Eq. (5.7): 

^5a s (m 2 )g(m 2 /Q 2 ) (5.10) 



5V 



o rrr 



where, setting m 2 /Q 2 = e, 



S(e) = -— Disc{.F(-e)} . 



(5.11) 



Since 5a s (m 2 ) is limited to low values of m 2 , the asymptotic behaviour of 5V at large 
Q 2 is controlled by the behaviour of jF(e) as e — > 0. We see from Eq. (5.11) that no 
terms analytic at e = can contribute to 5V. On the other hand for a square-root 
behaviour at small e, 



•F ~ a v — y/e 

271 



ay A\ 
dV = — , 

7T Q 



(5.12) 
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where Ai is the q — 1 moment of the strong coupling modification [13, 15], 

A, = ^ r^m«5a s (m 2 ). (5.13) 
Z7T Jo m 

Note that Sa s is a different quantity from the "effective coupling modification" 5a e s 
introduced in Ref. [16]. Their moments are related by the formula [15] 

C F f°° dm 2 2 sin(7rg/2) 

■m q 5a cS {m 2 ) = — v ' ' A q . (5.14) 



2tt Jo 



Thus in particular A-y = 2Ai/tt, and all the even moments of 5a e f[ are zero. 

We can express Ai in terms of the average value of a s in the infrared region, as 
follows [17, 18]. We substitute for 5a s in Eq. (5.13) 

5a s (m 2 ) ~ a s (m 2 ) — a^ T (m 2 ) , (5.15) 

where a" represents the expression for a s corresponding to the part already in- 
cluded in the perturbative prediction. As discussed in Ref. [17], if the perturbative 
calculation is carried out to second order in the MS renormalization scheme, with 
renormalization scale //|, then we have 

<*7(m 2 ) = «.04) + [bH^l/m 2 ) + k] a 2 M) (5-16) 
where for Nj active flavours (Ca = 3) 

b= ^ A -2N, (67-3^0,-HW, 

12tt ' 36tt v ; 

The constant k comes from a change of scheme from MS to the more physical scheme 
[19] in which a s is preferably defined at low scales. Then above some infrared match- 
ing scale /i : we assume that a s {m 2 ) and a"(m 2 ) approximately coincide, so that 

A± ~ — / dm (a s {m 2 ) - a s (nl) - [61n(//; \]m 2 ) + k\ a 2 (/4) ) 

TT JO ' 



where 



= (a^) - a s ^l) - [bHul/tf) + k + 2b] a 2 M)) , (5-18) 

cifo(/-0 = — / a s {m 2 ) dm . (5.19) 
Hi Jo 

Studies of event shapes in e + e~ annihilation [10, 11] suggest that «o — 0.5 for \i Y = 2 
GeV, which translates into a value of A x ~ 0.2 GeV for Q ~ 20 - 100 GeV. 

6. Power corrections to DIS event shapes 

We now use the method of the previous section to evaluate the coefficients of 1/Q 
for the DIS shape variables defined in Sect. 3. 
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6.1 Current jet mass 

Consider first the current jet mass, which is given in the soft region by the simplest 
expression, from Eq. (3.10): 

v(a,P)=aG(P-a) . (6.1) 
Eq. (4.9) therefore gives for the naive trigger function 

S%\e) = y/i, (6.2) 

dropping all terms of order e (i.e. order l/Q 2 ). 

In writing Eq. (6.1) we neglected the fact that the momentum fraction of the 
struck quark is not x but x(l + a). If we normalize to the Born cross section at fixed 
x, this means there should be an extra factor of q(x(l + a))/q(x) on the right-hand 
side, where q(x) is the relevant quark distribution. However, by Taylor expansion we 
see that this factor only gives a correction of order e, and so we may ignore it. 

Substituting Eq. (6.2) in Eq. (5.5), we obtain 

J#\e) = -2^ V~e ■ (6.3) 

71 

By simple dimensional analysis, the inclusive and non-inclusive corrections to the 
characteristic function, given by Eqs. (5.6) and (5.6), are also proportional to y/e: 

^\e) = -2^r mV Te (6.4) 



71 



where 



and 



where 



with 



°* r 2 ^ 2 (i+* 2 )] = 3 2994 c, 

A, Jo {i + x*){x + VTTx^) ' A)' 

^( e ) = -2— r raV /i (6.6) 

71 



m 



87T/3, 



,- / dk 2 t dz d<j) (C A M A + n f Mf) (k tl + k t2 - Jk 2 t +m 2 ) (6.7) 
o J 



k\ x = z[zk 2 + (1 — z)m 2 + 2\J z(l — z) k t m cos0] , 

k 2 t2 = (1 - z)[{\ - z)k 2 + zm 2 - 2y/z{l - z) k t m cos0] . (6.8) 

The integration in Eq. (6.7) was performed for e + e~ annihilation in Ref. [13]. As 
explained above, the result for DIS will be identical in the soft region, and therefore 
we can take the value directly from there (which we also confirmed numerically): 

rm = - J-(0.862CU + 0.052^) . (6.9) 

Po 
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Since the contributions to the jet mass from multiple soft parton emission are simply 
additive, see Eq. (3.9), there is no non-additive correction. The two-loop corrected 
characteristic function for the current jet mass is thus 

J^\e) = -2M— y/l (6.10) 

7T 

where M. is the Milan factor computed in Ref. [13], 

M = 1 + r in + r ni = 1 + (1.575CU - 0.104^)/ A, 

= 1.49 for n/ = 3 . (6.11) 

It follows from Eqs. (5.12) and (6.10) that the leading power correction to the 
mean current jet mass is 

4MAi A-, 
S (P) = —7V = L89 7T forn / = 3. (6.12) 

Comparing with our earlier result given in Ref. [1], 

<5<P>~2^, (6.13) 

we note the following differences: 

• The first moment of the effective coupling modification, Ai, is replaced by 
2Ai/tt where A\ is the first moment of the modification to the strong coupling 
itself, in accordance with Eq. (5.14). 

• The result is enhanced by a factor which is identical to the Milan factor obtained 
in e + e~ annihilation [13]. 

This means that universality of 1/Q corrections still holds, in the sense that the 
coefficient of 1/Q for the DIS current jet mass remains one half of that for the thrust 
in e + e~ annihilation, for example. Furthermore the numerical coefficient of A\/Q is 
only slightly lower than the coefficient of A\/Q calculated earlier using the "naive 
massive gluon" approximation, because the Milan factor M. is largely cancelled by 
the translation factor of 2 /it. Thus experimental results [9] interpreted earlier as 
measurements of Ai using the "naive" formula can be reinterpreted as giving a similar 
value for A\, taking into account the Milan and translation factors. 

6.2 Current jet thrust 

For the quantity r = 1 — Tq where Tq is the current jet thrust, we have from Eq. (3.3) 

v(a,P) = aG(P-a)+PG(a-P), (6.14) 
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and Eqs. (4.9) and (5.5) now give 

^\e) = 2V~e, rt\e) = -A^Ve. (6.15) 

7T 

Similarly, both the inclusive and non-inclusive corrections to the characteristic func- 
tion are twice as large as those for the jet mass, and again there is no non-additive 
part, so the leading power correction is 

= = 3 7Q Ai for 3 5 (6 16) 

n Q Q 

to be compared with our earlier result of AAi/Q in Ref. [1]. Once again the two-loop 
result (6.16) has a similar coefficient, but A\ is replaced hjA\. 

As stated earlier, the above result for the 1/Q power correction (unlike the per- 
turbative contribution) does not depend on whether the thrust is measured relative 
to the current axis or an axis which maximizes its value. We also see from Eq. (3.7) 
that the same leading power correction should be obtained if the thrust is normalized 
to the total energy in the current hemisphere, instead of Q/2. This is because the 
mean energy deficit in the current hemisphere, e, defined in Eq. (3.5), should have 
no 1/Q correction, since the two terms in Eq. (3.6) cancel on the average. 1 

6.3 C-parameter 

Similarly for the C-parameter (3.11), we have 

v(a,/?) = 12-^-6(0 -a) (6.17) 
a + p 

and the naive trigger function is 

nf\e) = 3ttv^, (6.18) 

giving 

jf ) (e) = -6C F y/e. (6.19) 

The inclusive and non-inclusive corrections scale in the same way and again there is 
no non- additive part, so 

5(C) = 12M ^ = 17.88^ for n f = 3 , (6.20) 

Our one-loop result for C in Ref. [1] was ambiguous, since it depended on how the 
gluon mass was included. Following the definition used here, the "naive" result is 
QttAi/Q = lS.SAi/Q. So again the two-loop result (6.20) has a similar coefficient, 
but with A\ replaced by A\. 



7. Conclusions 



x We expect the leading non-perturbative correction to (e) to be of order a s (Q 2 )/Q. 
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Our results on the coefficients of the leading power 
corrections to the DIS event shapes defined in Sect. 3 
are summarized in Table 1. For all the shape vari- 
ables studied, we have found that, as in e + e~ anni- 
hilation [13], the enhancement factor for the leading 



V 


1-loop 2-loop 


(p) 


2 1.9 


(r) 


4 3.8 


(C) 


6tt - 12tt 17.9 



power correction is equal to the Milan factor (6.11), Table 1: Coefficients of A x jQ 
provided the "naive" coefficient is computed in the (Hoop) and A±/Q (2-loop). 
way specified in Sect. 4. This eliminates the ambigu- 
ity in the prediction for the C-parameter, which was present at the 1-loop level [1] 
as indicated. 

In all cases the enhancement factor practically cancels the suppression factor of 
2/ Ti that comes from translating the first moment of the effective coupling modifi- 
cation, Ai, into the corresponding quantity for the strong coupling itself, A\. This 
quantity can be used to infer the value of do, the low-energy moment of a s defined 
in Eq. (5.19). 

The experimental data on thrust and jet masses in both DIS [9] 2 and e + e~ 
annihilation [10, 11] suggest a value of do ~ 0.5. Final experimental results on the 
C-parameter in DIS have yet to be presented. 
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